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Method and system for measuring in a dynamic sequence of medical Images 

Description 

5 Technical field 

The present invention relates to a method and a system for measuring in a 
dynamic sequence of medical images of a moving body part. 

10 Background pf iilJB invention 

Non-Invasive methods for measuring functions and dysfunctions of body 
organs typically involve scanning with an X-ray device to obtain projection 
data for generation of a series of images. When measurements are to be 

15 performed of a moving body part, Images of several events of the movement 

have to be created. Images are then generated with a certain frequency, 
whereby a sequence of consecutive images is formed. The sequence depicts 
the same anatomical area over a period of time during which some dynamic 
event takes place and is thus a dynamic sequence. An example of such a 

20 dynamic sequence is an image sequence of the heart over one or more 

heartbeat cycles. 

In order to be able to perform measurements in the images, different kinds 
of measurement tools for application to the images are needed- A 
25 fundamental tool in image measurements is the distance measurement tool. 

Another important tool Is the angle measurement tool for measuring the 
angle between two lines superimposed on an image. 

Common for many tools available today is that they measure static data in 
30 static images. When a distance measurement is to be done in a static image, 

a user identifies two anatomical features in the image and applies a distance 
measurement tool to determine the distance between them. 

There are tools and methods for measuring In dynamic sequences, but they 
: 35 consist of two or more instances of an intrinsically static tool that is more or 

: I : less manually applied to several of the Images in a dynamic sequence. The 

result from these different instances can then be used for comparisons or to 
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compute cruder approximations to dynamic entities by using ratios or Huvudfoxen Ka 
differentiations. 

For example, a commonly used method today for measuring in a dynamic 
5 sequence is distance measurement applied to an image sequence from 

coronary angiography Involving the following steps: manual Identification of 
two images in the sequence, for example showing the end-diastole and end- 
systole phase In the same heart cycle; manual location of the most apical 
part of the coronary artery branches in both images; manual location of the 

10 lower contour of the left coronary ostium In both images; manual location of 

two points on the horizontal part of the circumflex artery, one proximal and 
one distal, in both images; measurement of the distance from the point from 
the second step to the point from the third step In one image and comparison 
with the distance between the same points In the other Image; and 

15 measurement of the distance from the point from the second step to the 

points from the fourth step and comparison with the distance between the 
same points in the other image. 

Performing the above mentioned steps only once gives very limited dynamic 
20 information since only two images in the sequence are used for the 

measurements. To obtain more dynamic information, a user has to repeat 
the steps for every image in the dynamic sequence or at least many images 
In the dynamic sequence. Usually there are over 100 images In a sequence. 
Thus, using the above mentioned method on all images in an image 
25 sequence is rather time-consuming. 

Through EP 1 088 517 a method and an apparatus for motion-free cardiac 
imaging is known. In EP 1 088 517 a fixed reference point and a dynamic 
point are used for selecting images without motion-induced artefacts in a 
30 dynamic sequence of computerized tomographic images, A line is drawn from 

_ : the reference point to the dynamic point in all of the Images. The line length 

represents the distance between the reference point and the dynamic point 
: and an Image or images are selected in which the line length remains 

constant relative to the previous image. 
: 35 Image recognition software can be used either to identify a pair of reference 

: points for measurements of relative motion or, once reference points are first 

identified, to identify corresponding reference points on other Images. 
Distances between the automatically identified points can be selected by 
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software, based upon selected criteria. However/ in EP 1 088 517 Huvudfaxeo Kosscn 
determination of other variables such as speed, acceleration and retardation 
are not described for measurement in a dynamic sequence of images of a 
moving body organ. Furthermore, only two points, one fixed reference point 
and one dynamic point are used, it is an advantage to be able to use more 
than two points as well as more than one dynamic point. For example, angle 
measurements and area measurements require at least three points and 
measurements of relative motion of two moving parts of a moving body part 
require at least two dynamic points. 



Accordingly, it is an object of the present Invention to provide an improved 
method for measuring In a dynamic sequence of medical images of a moving 
body part. 

This object is achieved by means of a method comprising the steps of: 
scanning a portion of a body of a patient including the moving body part for 
generating time resolved projection data; generating said images from said 
projection data; defining at least one measurement point in the moving body 
part in one of said images; defining a reference point In one of said images to 
a point being fixed relative to the Image geometry, said reference point being 
different from said at least one measurement point; automatically tracking 
the at least one measurement point in all of said images of the sequence; 
automatically indicating the reference point in ail of said images of the 
sequence; automatically determining a length and a direction of a vector 
extending from the reference point to one of the at least one measurement 
points for each pair of reference point and one measurement point in ail of 
said images of the sequence and, automatically determining a rate of change 
of said length and/or said direction of said vector{s) between selected images 
In said sequence of images. 

Another object of the present invention is to provide an Improved system for 
measuring in a dynamic sequence of medical images of a moving body part. 



This object Is achieved through a system for measuring In a dynamic 
sequence of medical images of a moving body part, said system having 
means for: scanning a portion of a body of a patient including the moving 
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body part for generating time- resolved projection data; generating said 
images from said projection data; defining at least one measurement point in 
the moving body part in one of said Images; defining a reference point In one 
of said images to a point being fixed relative to the image geometry, said 
5 reference point being different from said at least one measurement point; 

automatically tracking the at least one measurement point in all of said 
images of the sequence; automatically Indicating the reference point in all of 
said Images in the sequence; automatically determining a length and a 
direction of a vector extending from the reference point to one of the at least 
10 one measurement points for each pair of reference point and one 

measurement point in all of said images of the sequence and, automatically 
determining a rate of change of said length and/or said direction of said 
vector(s) between selected images in said sequence of images. 

15 Advantages of a method and a system according to the present invention 

include possibility to automatically measure movement variables such as 
speed, acceleration/retardation and/or direction of movement of a point or 
points of the moving body part, to use more than two points for 
measurement, to use more than one dynamic point for measurement, to 

20 automatically measure angles and areas and to use one-dlmenslonal, two- 

dimensional and/or time resolved two-dimensional search fields for the 
automatically tracking of the dynamic pomt(s). 

Still other objects and advantages of the present invention will become 
25 apparent from the following detailed description considered In conjunction 

with the accompanying drawings. It is to be understood, however, that the 
drawings are designed solely for purposes of Illustration and not as a 
definition of the limits of the invention, for which reference should be made 
to the appended claims. It should further be understood that the drawings 
30 are not necessarily drawn to scale and that, unless otherwise indicated, they 

are merely Intended to conceptually illustrate the methods and systems 
herein. 



35 
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In the following, the Invention Is described In greater detail applied by way of 
example to a heart. The description is made with reference to attached 
drawings, in which like reference characters denote similar elements and 

Fig. l shows a composite of Images In a dynamic sequence of schematic 
images of parts of a heart with a fixed reference point set in an apical part of 
the heart and a dynamic point set In a moving part of the left ventricular 
wall; 

Fig. 2a-d show separate images of the dynamic sequence shown In fig. 1 and 
are also views of distance measurement; 

Fig. 3a-d show a first embodiment of the invention including the images 
shown in fig. 2a-d; 

Fig. 4a-d show a second embodiment of the invention including the images 
shown in fig. 2a-d; 

Fig. 5a-d show the images shown in fig. 2a-d, but with a second dynamic 
point set in the left ventricular wall; 

Figure 6a-d are views of measurement of dynamic angles and dynamic areas 
In the images shown In figure 5a-d. 

Fig. 7a-d show a fourth embodiment of the present invention including the 
Images shown in ftg. 5a-d; 

Fig. 8a-d show a fifth embodiment of the present invention including the 
images shown in fig. 5a~d; 

Fig. 9a-d show a sixth embodiment of the present invention Including the 
images shown in fig. 5a-d; 



Figure 10 shows a schematic view of four points used for angiographic 
measurements. 
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5 Figure l shows a composite of highly schematic Images of some parts of a 

heart 1. The images shown In figure 1 belong to a dynamic sequence of 
images of a heart 1, but are only a few images of a dynamic sequence 
covering a cardiac cycle of the heart 1. The sequence of images is generated 
from projection data from time resolved two dimensional X-ray scanning of a 
10 portion of a body of a patient including the heart 1. The projection data Is 

thus time resolved two-dimensional data. In order to cover a complete 
cardiac cycle the images are generated at an adequate frequency , for 
example 12/5 Images per second. A contrast medium is introduced before the 
scanning to visualize vessels properly on the images. 

15 

The heart 1 is a cyclically moving body part having a right atrium 2, a left 
atrium 3, a right ventricle 4 and a left ventricle 5. A single cardiac cycle 
consists of one diastolic phase (expansion phase) and one systolic phase 
(contraction phase) for the atria 2,3 and one diastolic phase and one systolic 

20 phase for the ventricles 4,5, The right ventricle 4 receives deoxygenated 

blood during its diastolic phase and pumps it into the pulmonary arteries 
during Its systolic phase, while the left ventricle 5 receives oxygenated blood 
during its diastolic phase and pumps It Into the aorta during Its systolic 
phase. The heart muscles are supplied with blood by the coronary arteries, 

25 which encircle the heart and thus follow the expansion and contraction of the 

heart 1 during a cardiac cycle. 



For different reasons it would be advantageous to be able to track a point of 
the heart 1 and measure for example movement variables of that point of the 
30 heart 1 relative a fixed point during the cardiac cycle. Movement variables, 

which might be Interested to measure are distance, speed, acceleration, 
retardation, direction of movement of a point of the heart as well as dynamic 
angles and areas. 

35 Using the method according to the invention, a reference point 6 and a 

measurement point 7 have to be defined by an user. The reference point 6 Is 
set to a point which is fixed in its position relative the image geometry, \ e it 
is set to a point of a structure which is not essentially moved during a cardiac 
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cycle. In the images shown in figure 1 the reference point 6 is set to a n Huvu ^ 2X2n Kossan 
apical part of the heart l, which is not essentially moved during a cardiac 
cycle. The reference point 6 is set in one Image by a user and Is then 
automatically Indicated, i e marked, in the other images in the sequence 
5 through image processing software. 

In the Images in figure 1 a point of the wall of the left ventricle 5 is desired 
to be measured and then that point Is defined as the measurement point 7. 
The measurement point 7 is a dynamic point, since it Is moved during the 
cardiac cycle. The dynamic point 7 Is also set by an user in one Image, and 

10 the corresponding points In the other Images are then automatically tracked 

by using image processing software including one of a number of known 
algorithms for automatic tracking of anatomic parts. An example of one such 
algorithm Is FMI-SPOMF (Fourier-Mellin Invariant Symmetric Phase-Only 
Matched Filtering), described for instance by Qin-sheng Chen, Michel Defrise 

15 and F. Deconick in IEEE Transactions on Pattern Analysis and Machine 

Intelligence, Vol. 16, No. 12, pp 1146-1168, Dec 1994. In order to improve 
the efficiency of the algorithm/algorithms used, the images can be pretreated 
to increase the contrasts between the imaged object and the background. 
The tracking Is performed using information about the imaged object and the 

20 darkest point can for example be tracked since It probably represents a 

vessel containing a lot of contrast medium. 



The wall of the left ventride 5 expands during the diastolic phase and 
contracts during the systolic phase. Thus, during the diastolic phase of the 
25 left ventricle 5 the dynamic point 7 is moved from a position denoted 8 Into 

expanded positions denoted 9, 10 and 11. The position 8 represents the end- 
systolic position of the point 7 and the position 11 represents the end- 
diastollc position of the point 7. Accordingly, during the systolic phase of the 
left ventricle 5, the point 7 is moved from the position 11 into position 8. 

30 

Figures 2a-d show separate images of the images in the dynamic sequence 
shown in fig. i and are also views of distance measurement. The Images In 
figures 2a-d are subsequent In time. Figure 2a shows an end-systolic view of 
the heart 1 and figure 2d an end-dlastolic view of the heart 1. A length and a 
35 direction of a vector v a extending from the reference point 6 to the first 

dynamic point 7 can be automatically determined in all of the images of the 
sequence using Image processing software. A first distance d n between the 
dynamic point 7 and the reference point 6 can be also be automatically 
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determined In each of the images In the sequence of images using the 
of the vector v t , The first distance dn between the reference point 6 and the 
dynamic point 7 can also be automatically compared between different 
images and used to determine the change of the position, i e the movement, 
5 of the dynamic point 7 during the cardiac cycle. 

Since the data used for reconstruction of the Images in figure 1 is time 
resolved two-dimensional data, it is also possible to determine a rate of 
change of the first distance dn, I e movement variables, between different 
10 images in the sequence of image. The rate of change is also automatically 

determined using image processing software. 

Movement variables which can be determined are for example speed, 
acceleration, retardation and direction of movement of the dynamic point in 
different moments of the cardiac cycle, i e In different images of the 
15 sequence of images. Furthermore, the speed, acceleration and/or retardation 

with which the dynamic point 7 is moved between its different positions In 
two different images can also be determined. The movement variables can be 
measured in all of the below described embodiments of the present 
Invention. 

20 

Figures 3a-d show a first embodiment of the invention including the images 
shown In figures 2a-d. In the first embodiment of the invention a one- 
dimensional search field is used for the tracking of the dynamic point 7 in the 
different images In the sequence of Images. A vector 12, i e a line, 
25 representing the one-dimensional search field is marked In one of the images; 

and the dynamic point 7 is then automatically tracked along that line 12 in all 
of the images in the sequence of Images using Image processing software as 
previously mentioned. 



30 Figures 4a-d show a second embodiment of the invention including the 

images shown In figures 2a-d. In the second embodiment of the invention a 
two-dimensional search field is used for the tracking of the dynamic point 7 
in the different images in the sequence of images. An area 13, i e a 
rectangle, representing the two-dimensional search field Is marked In one of 

35 the images and the dynamic point 7 Is then automatically tracked within that 

rectangle 13 in all of the images in the sequence of images using image 
processing software as previously mentioned. 
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In a third embodiment (not shown) of the invention a time resolved two- 
dimensional search field Is used for the tracking of the dynamic point 7 In the 
different Images In the sequence of images. An area 13, I e a rectangle, 
representing the two-dimensional search field is marked in one of the images 
5 In the sequence of Images and the dynamic point 7 is automatically tracked 

within the rectangle using image processing software as previously 
mentioned, but information from previous and subsequent images are used 
as guide to which point of possible points in the rectangle that is the 
searched dynamic point 7. Expected positions of the dynamic point 7 based 
10 on information from previous images can also be used to choose which point 

in the rectangle that Is the searched dynamic point 7. 



Figures 5a-d show the images shown in figures 2a-d, but with a second 
dynamic point 14. Using the method according to the invention/ movement of 

15 two different parts of the moving body part can simultaneously be measured 

by defining a second dynamic point 14, which is different from the first 
dynamic point 7. In the figures 5a-d the second dynamic point 14 Is set to a 
point of the wall of the left ventricle 5 being different from the first dynamic 
point 7. The second dynamic point is moved from a position denoted 19 into 

20 expanded positions 20, 21 and 22. The position 20 represents the end- 

systolic position of the point 14 and the position 22 represents the end- 
diastoiic position of the point 14. Accordingly, during the systolic phase pf the 
left ventricle 5, the point 14 is moved from the position 22 into position 19. A 
length and a direction of the vector v x extending from the reference point 6 

25 to the first dynamic point 7 as well as a vector v 2 extending from the 

reference point 6 to the second dynamic point 14 can then be automatically 
determined using image processing software. 

The first distance d ±1 between the first dynamic point 7 and the reference 
point 6 as well as the first distance 6 12 between the second dynamic point 14 

30 and the reference point 6 can then also be automatically determined using 

the length of the vector v t and the vector v 2 respectively and compared 
between different images. Furthermore, the movement between the position 
in one image and the position In another image of both the first dynamic 
point 7 and the second dynamic point 14 as well as movement variables like 

35 for example the above mentioned can also be determined. It is also possible 

to automatically determine the movement of the two dynamic points 7,14 
relative each other through determining a second distance d 2 i between the 
two dynamic points 7,14 and comparing the second distance d 2 i between 
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different images. The second distance d 2 x Is automatically determined using 
the lengths of the vectors v x and v 2 and using image processing software. 
Furthermore, using two dynamic points 7,14 it Is also possible to measure 
dynamic angles and areas. Dynamic angles and areas can be measured in all 
of the below described embodiments. 

Figures 6a-d are views of measurement of dynamic angles and dynamic 
areas in the images shown in figures 5a-d. A dynamic angle x can be 
automatically determined by using image processing software and by using 
the first distances dn, d i2 and the second distance d 2 iin each of the images 
of the sequence. Other angles than the shown angle x can of course also be 
determined. A dynamic area a can automatically be determined by using 
image processing software and by using the first distances du, d i2 and the 
second distance d 2 i in each of the images of the sequence. 

Figures 7a-d show a fourth embodiment of the Invention including the images 
shown in figures 5a-d. In the fourth embodiment of the Invention one- 
dimensional search fields are used for the tracking of the first and second 
dynamic points 7,14 respectively in the different Images in the sequence of 
images. A vector 12, i e a line, representing the one-dimensional search field 
is marked for each dynamic point 7,14 in one of the images and the dynamic 
points 7,14 are then respectively tracked along the lines 12 in all of the 
images in the sequence of images using images processing software as 
previously mentioned. 

Figures 8a-d show a fifth embodiment of the invention including the images 
shown in figures 5a-d. In the fifth embodiment of the invention two- 
dimensional search fields are used for the tracking of the first and second 
dynamic points 7, 14 respectively In the different Images in the sequence of 
images. An area 13, 1 e a rectangle, representing the two-dimensional search 
field is marked for each dynamic point 7,14 in one of the images and the 
dynamic points 7 f 14 are then respectively tracked within that rectangles 13 
In all of the Images in the sequence of images using Image processing 
software as previously mentioned. 



Figures 9a-d show a sixth embodiment of the invention including the images 
shown in figures 5a-d. In the sixth embodiment of the Invention a one- 
dimensional search field is used for the tracking of one of the first and the 
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second dynamic points 7,14 and a two-dimensional search field Is usec^for ^°' <en Kelson 
the tracking or the other of the first and the second dynamic points 7,14 in 
the different Images In the sequence of Images. A vector 12, I e a line 
representing the one-dimensional search field and an area 13, 1 e a 
5 rectangle, representing the two-dimensional search field are marked In one 

of the images and the dynamic points 7,14 are then tracked along the line 
and within the rectangle respectively In all of the images in the sequence of 
Images using Images processing software as previously mentioned. In the 
figures 9 a-d a one-dimensional search field is used for the tracking of the 
10 first dynamic point 7 and a two-dimensional search field Is used for the 

tracking of the second dynamic point 14, but Instead a one-dimensional 
search field can be used for the tracking of the second dynamic point 14 and 
a two-dimensional search field can be used for the tracking of the first 
dynamic point 7. 

15 

In a seventh embodiment (not shown) of the invention time resolved two- 
dimensional search fields are used for the tracking of the dynamic points 
7,14 respectively in the different images in the sequence of images. An area 
13, i e a rectangle, representing the two-dimensional search field Is marked 

20 for each dynamic point 7,14 In one of the images in the sequence of images 

and the dynamic points 7,14 are respectively tracked within the rectangles 
using images processing software as previously mentioned, but information 
from previous and subsequent images are used as guide to which points of 
possible points In the rectangles that are the searched dynamic points 7,14. 

25 Expected positions of the dynamic points 7,14 based on information from 

previous images can aiso be used to choose which points are the searched 
dynamic points 7,14. 



In an eighth embodiment (not shown) of the Invention a one-dimensional 
30 search field is used for the tracking of one of the first and the second 

dynamic points 7,14 and a time resolved two-dimensional search field is used 
for the tracking of the other of the first and the second dynamic point 7,14 in 
the different images in the sequence of Images. A vector 12, I e a line 
representing the one-dimensionai search field and an area 13, i e a 
35 rectangle, representing the two-dimensional search field are marked in one 

of the images and the dynamic points 7,14 are then respectively tracked 
along the line 12 and within the rectangle 13 in all of the Images In the 
sequence of images using image processing software as previously 
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mentioned. However, information from previous and subsequent images are 
used as guide to which point of possible points in the rectangle 13 that Is the 
searched dynamic point 7,14. Expected positions of the dynamic point 7,14 
based on information from previous images can also be used to choose which 
5 point in the rectangle Is the searched dynamic point 7,14. 

In a ninth embodiment (not shown) of the invention a two-dimensional 
search field is used for the tracking of one of the first and the second 
dynamic points 7,14 and a time resolved two-dimensional search field Is used 

10 for the tracking of the other of the first and the second dynamic points 7,14 

In the different images in the sequence of images. An area 13, i e a 
rectangle, representing the two-dimensional search field is marked for each 
dynamic point 7,14 in one of the images and the dynamic points 7,14 are 
then respectively tracked within that rectangles 13 in all of the images in the 

15 sequence of images using image processing software as previously 

mentioned. However, information from previous and subsequent images are 
used as guide to which point of possible points in the rectangle 13 
representing the time resolved two-dimensional search field that is the 
searched dynamic point 7,14. Expected positions of the dynamic point 7,14 

20 based on information from previous Images can also be used to choose which 

point in the rectangle 13 representing the time resolved two-dimensional 
search field that Is the searched dynamic point 7,14. 

It must be obvious that the method according to the Invention can be 
25 expanded to cover more than two dynamic points as well as more than one 

reference point. Any combinations of one-dimensional, two-dimensional and 
time resolved two-dimensional search fields for the tracking of the different 
points are then possible to use. 

. 30 By way of example the method according to the invention is described 

involving coronary angiography. Projection data of the coronary arteries are 
generated through coronary angiography, i e X-ray examination after 
injection of an opaque dye. The Images are generated during at least one 
cardiac cycle at a frequency of 12.5 Images per second. The right anterior 
35 oblique 30° projection has shown good results and is preferably used. A 

sequence of images are then generated through reconstruction from the 
projection data. Four points 15, 16, 17, 18 are then located and marked In 
one image in the sequence of Images; a point 15 representing the most 
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apical point of the coronary artery branches, a point 16 representing the 
lower contour of the left coronary ostium, a proximal point 17 and a distal 
point 18 on the horizontal part of the circumflex artery. The corresponding 
four points 15, 16, 17, 18 are then automatically located In each of the 
5 Images In the image sequence using image processing software. Figure 10 

shows a schematic view of the four points 15, 16, 17, 18 used for 
angiographic measurements. 

The eplcardlal part of the apex is almost stationary during the heart cycle. 
10 Thus, the most apical point 15 of the coronary branches is almost stationary 

during the heart cycle and is therefore used as a reference point. The other 
three previously mentioned points 16, 17, 18 are dynamic points 
representing a first, a second and a third dynamic point. A one-dimenslonal, 
a two-dimensional or a time resolved two-dimensional search field can be 
15 used for the tracking of each of the corresponding dynamic points 16, 17, 18 

in the other Images In the sequence of Images. The first distance du from the 
most apical point 15 of the coronary branches to the point 16 representing 
the lower contour of the left coronary ostium Is then automatically 
determined In each Image showing the amplitude of the left coronary ostium 
20 motion. The first distance d i2 between the most apical point 15 of the 

coronary branches and the proximal point 17 on the horizontal part of the 
circumflex artery as well as the first distance d*j between the most apical 
point 15 of the coronary branches and the distal point 18 on the horizontal 
part of the circumflex artery are also automatically determined in each Image 
25 showing the motion of the proximal point 17 and the distal point 18 

respectively. Furthermore, the second distance dm between the proximal 
point 17 and the distal point 18 of the circumflex artery is determined, which 
second distance d 22 3 represents a portion of the circumflex artery extending 
along the atrioventricular groove. The motion of that portion of the 
: 30 circumflex artery roughly follows the most basal part of the left ventricular 

* . wail, which means that change of the distance between the proximal and the 

distal point of the circumflex artery represents the shortening and 
: lengthening of the left ventricle 5. 

35 Although cardiac Imaging embodiments are described in detail, the Invention 

has more general applicability and can be applied to another moving body 
part, such as for example the lungs. Furthermore, the method according to 
the invention can also be expanded to work on time resolved three- 
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dimensional angiographies as well as simultaneous measuring of dynamic 
points in different body parts. 

Thus, while there have been shown and described and pointed out 
5 fundamental novel features of the invention as applied to embodiments 

thereof, it will be understood that various omissions and substitutions and 
changes In details of the methods described, and In their operation, may be 
made by those skilled in the art without departing from the spirit of the 
invention. For example/ it Is expressly intended that all combinations of those 

10 method steps and/or system elements which perform substantially the same 

function In substantially the same way to achieve the same results are within 
the scope of the invention. Moreover, it should be recognized that method 
steps and/or system elements shown and/or described in connection with 
any disclosed form or embodiment of the invention may be incorporated in 

15 any other disclosed or described or suggested form or embodiment as a 

general matter of design choice. It is the intention, therefore, to be limited 
only as indicated by the scope of the claims appended hereto. 





1 
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1. A method for measuring in a dynamic sequence of medical images of a 
moving body part comprising the steps of: 

5 - scanning a portion of a body of a patient including the moving 

body part for generating time resolved projection data; 
generating said images from said projection data; 
defining at least one measurement point (7,14) in the moving 
body part in one of said images; 
10 - defining a reference point (6) in one of said Images to a point 

being fixed relative to the image geometry, said reference point 
(6) being different from said at least one measurement point 
(7,14); 

- automatically tracking the at least one measurement point (7,14) 
15 in all of said images of the sequence; 

automatically indicating the reference point (6) in all of said 
images of the sequence; 

automatically determining a length and a direction of a vector (v i/ 
v 2 ) extending from the reference point (6) to one of the at least 
20 one measurement points (7,14) for each pair of reference point 

(6) and one measurement point (7,14) in all of said images of the 
sequence and, 

automatically determining a rate of change of said length and/or 
said direction of said vector(s) (vi, v 2 ) between selected Images In 
25 said sequence of images. 

2. A method according to claim l, which method further comprises the step 
of automatically determining a first distance (d U/ d 12 ) between the 
reference point (6) and one of the at least one measurement points 

30 (7,14) for each pair of one reference point (6) and one measurement 

point (7,14) using the length of the corresponding first vector {v u v 2 ). 

3. A method according to any of the preceding claims, which method further 
comprises a step of automatically determining a direction of movement of 

35 the at least one measurement point (7,14) by using the direction of the 

corresponding first vector (y x , v 2 ). 
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4. A method according to any of the preceding claims, which method further 
comprises a step of automatically determining a speed of the at least one 
measurement point (7,14) by using said rate of change of the length of 
the corresponding first vector (vj, v 2 ). 

5 

5. A method according to any of the preceding claims, which method further 
comprises a step of automatically determining an acceleration and/or a 
retardation of the at least one measurement point (7, 14) by using said 
rate of change of the length of the corresponding first vector (v t , v 2 ). 

10 

6. A method according to any of the preceding claims, wherein the step of 
defining at least one measurement point (7,14) In the moving body part 
in one of said Images comprises defining at least two measurement points 
(7,14). 

15 

7. A method according to claim 6 f wherein the method further comprises the 
step of automatically determining a second distance (d 2 0 between two of 
the at least two measurement points (7,14) for each pair of two 
measurement points (7,14) using the first lengths of the corresponding 

20 first vectors (v 1# v 2 ). 



8. A method according to any of the preceding claims, which method further 
comprises a step of automatically comparing said rate of change of said 
length and/or said direction of said vector(s), said first distance(s) (dn, 

25 d 12 ), said speed(s), said acceleration(s)/reterdation(s), said direction of 

movement(s) and/or said second distance(s) (d 21 ) between selected 
images in said sequence of images. 

9. A method according to any of the preceding claims, wherein the step of 
30 automatically tracking the at least one measurement point (7,14) is 

preceded by a step of creating in each of said images a one-dimensional 
search field for the tracking of each of the at least one measurement 
points (7,14). 



35 



10. A method according to any of claims 1-8, wherein the step of 

automatically tracking the at least one measurement point (7,14) is 
preceded by a step of creating in each of said images a two-dimensional 
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search field for the tracking of each of the at least one measuremiw J ' 

polnts(7,14). 

11. A method according to any of claims 1-8, wherein the step of 
5 automatically tracking the at least one measurement polnt(7,14) is 

preceded by a step of creating In each of said images a time resolved 
two-dimensional search field for the tracking of each of the at least one 
measurement points (7,14). 

10 12. A method according to claim 11, wherein the step of creating In each of 

said images a time resolved two-dimensional search field for the tracking 
of the at least one measurement point (7 r 14) further comprises creating 
the search field using information from previous and/or following Images 
in said sequence of Images. 

15 

13. A method according to claim 11, wherein the step of creating in each of 
said images a time resolved two-dimensional search field for the tracking 
of the at least one measurement point (7,14) further comprises creating 
the search field using expected values based on Information from the 
20 previous Image in said sequence of images. 



14. A method according to any of claims 7-13, which method further 
comprises a step of determining at least one dynamic angle using said 
first distance(s) (dn, di2) and/or said second distance(s) (d zl )- 

25 

15. A method according to any of claims 7-14, which method further 
comprises a step of determining an area using said first distances (dn, 
di2) and/or said second distances (d 21 ). 

30 16. A method according to any of claims 7-15, wherein the step of 

automatically tracking the at least two measurement points (7,14) is 
preceded by a step of creating in each of said Images any combinations of 
one-dimensional search fields, two-dimensional search fields and time 
resolved two-dimensional search fields for the tracking of the at least two 

35 measurement points (7,14)» 



17. A method according to any of the preceding claims, wherein the step of 
automatically tracking the measurement point(s) (7,14) Is preceded by a 
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step of Image processing increasing the contrast between the reproduced 
object(s) and the background. 

18, A method according to any of the preceding claims, wherein the step of 

5 automatically tracking the measurement point(s) (7,14) comprises image 

processing software comprising at least one algorithm tracking the 
measurement polnt(s) 

19, A system for measuring in a dynamic sequence of medical images of a 
10 moving body part, said system having means for: 

scanning a portion of a body of a patient including the moving 
body part for generating time resolved projection data; 
generating said images from said projection data; 
defining at least one measurement point (7,14) in the moving 
15 body part in one of said images; 

- defining a reference point (6) in one of said Images to a point 
being fixed relative to the image geometry, said reference point 
(6) being different from said at least one measurement point 
(7,14); 

20 - automatically tracking the at least one measurement point (7,14) 

in all of said images of the sequence; 

automatically indicating the reference point (6) In all of said 
images of the sequence; 

automatically determining a length and a direction of a vector (v lr 
25 v 2 ) extending from the reference point (6) to one of the at least 

one measurement points (7,14) for each pair of reference point 
(6) and one measurement point (7,14) in ail of said images of the 
sequence and, 

- automatically determining a rate of change of said length and/or 

: 30 said direction of said vector(s) (v l7 v 2 ) between selected images in 

said sequence of Images. 

20. A system according to claim 19, wherein the scanning is performed by 
a X-ray device. 

35 
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Abstract 



The present invention relates to a method and a system for measuring in a 
dynamic sequence of medical Images of a moving body part. A reference point 
5 being fixed relative to the image geometry and at least one measurement point 
are defined in the moving body part in one Image In the sequence of images. 
The reference point Is then automatically indicated and the at least one 
measurement point is then automatically tracked in ail of said images of the 
sequence. A length and a direction of a vector extending from the reference 
10 point to one of the at least one measurement points for each pair of reference 
point and one measurement point is automatically determined In all of said 
images of the sequence and a rate of change of said length and/or said direction 
of said vector(s) is automatically determined between selected images in said 
sequence of images. 



15 
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Fig. 2a 
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